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Abstract. Recent OGLE (Optical Gravitational Lensing Experiment) and GLITP (Gravitational Lens 
International Time Project) monitoring data for QSO 2237+0305 (Huchra et al. 119851 have been analyzed through 
a newly optimized N-body microlensing analysis method, the Local HAE Caustic Modeling (LOHCAM). This 
method simultaneously solves for the size of the source and N-body HAE (High Amplification Events) caustic 
shapes in the source plane and determines those sizes only as a function of the projected transverse velocity of the 
source. By applying this method to the light curves of the A & G lensed components in the Einstein Gross, these 
data are accurately reconstructed for the first time. From these modeling studies, we report several interesting 
results: the minimum number of microlenses required for possible caustic models, the possible evidence of an 
accretion disk with a central hole at the heart of the quasar, the size of the UV-continuum source, the masses of 
the microlenses being directly responsible for the observed HAEs, the estimated mass range of a super massive 
black hole (SMBH) in QSO 2237+0305 and finally some clues for the direction of the source motion in the sky. 
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1. Introduction 

Gravitational microlensing effects induced by compact ob¬ 
jects in a galaxy can be used as a cosmic microscope 
to probe at micro-arcsec angular scales the structure of 
the lensed background quasar (Chang & Refsdal II979I 
Young II98II Chang & Refsdal II984L Kayser, Refsdal & 
Stabell [T98till . A few years ago, the OGLE team (Wozniak 
et al. iinnnj has reported impressive microlensing light 
curves for the 4-macro images (hereafter denoted A, B, 
C & D) of QSO 2237+0305 {z ~ 1.695) lensed by a 
nearby spiral galaxy {z ^ 0.0394). This system consti¬ 
tutes a very promising candidate for microlensing studies 
(Kayser, Refsdal & Stabell lTOHOL Kayser & Refsdal fl989ll . 
Indeed, because the expected time delays are so short, 
about one day (Schneider et al. 1 198511 . intrinsic variability 
should show up almost simultaneously in the four images 
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SO that any brightness variations affecting only some of 
the four single images are necessarily due to microlensing. 
Also, because of the large distance ratio between the QSO 
and the lens, microlensing should lead to more frequent 
and rapid high amplification events (HAEs) than for any 
other known gravitational lens system. So far, there have 
been several monitoring observations for this unique grav¬ 
itational lens system, however, the OGLE data (Wozniak 
et al. l2()()()l) are unique in that they cover prominent HAEs 
with unprecedented daily time resolution. Recent GLITP- 
ISIS data (Alcalde et al. I2()()2[ Moreau et al. 1200411 pro¬ 
vide some additional data around the microlensing peak 
of A and C. In section 2, we present data handling for the 
OGLE and GLITP data. In section 3, the LOHCAM tech¬ 
nique is introduced. In section 4, we present the LOHCAM 
results and discussion. 
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Fig. 1. OGLE+GLITP data and model fitting results. Combined V-band light curves for image A and the possible 
caustic model, (a) The OGLE (231 points) & GLITP-ISIS (53 points; around day 1500) data are shown altogether 
with the simulated light curve (solid line) using a disk with a central hole source model. Note the smooth and fine 
htting over the entire data points, especially around the HAE peak (see more details in Fig. 3). (b) The corresponding 
HAE caustic network (solid line) is shown with the straight track (dotted lines) for the moving source (Y=0). The 
concentric circles represent the sizes of the best-htted inner and outer radii of the face on accretion disk. The sizes of 
the caustic network and of the source are represented in time units (days) along the X and Y axes. 


2. Data handling 

The OGLE data cover approximately 3 years with 231 
data points, and the GLITP data cover 120 days with 
53 data points. We note that two different methods have 
been applied to GLITP data analysis: GLITP-PSF fitting 
data & GLITP-ISIS data (Alcalde et al. 120021 Moreau 
et al. Eoni. In this study, we use the OGLE data and 
GLITP-ISIS data which provided smaller error bars com¬ 
pared to GLITP-PSF fitting data. In order to combine 
these two data sets, first we corrected a zeropoint differ¬ 
ence of 0.065 and 0.13 mag. in the V band between the 
OGLE and GLITP-ISIS data of A & G, respectively, so as 
to match the slope of the OGLE data. Then, the original 
logarithmic magnitude scale is transformed into a linear 
flux unit (mjy) based on the same formulation given by 
the OGLE team. Figures la and 2a illustrate the resulting 
combined light curves in mJy unit that we have analyzed. 

3. Local HAE Caustic Model 

Usefulness of lightcurves of quasar microlensing is studied 
and discussed by many authors (Wambsganss Paczynski 
& Schneider I1990L Agol & Krolik [T999I and Mineshige & 
Yonehara If 009|l . They have investigated various methods 
to reconstruct spatial structure of quasar accretion disk. 
Recently, several authors fYonehara l2()()ll ShalvaDin l2()()ll 
Shalyapin et al. Enni have attempted to fit parts of the 


OGLE or GLITP data using HAE approximation formu¬ 
lae for a fold caustic and a cusp under the assumption 
of a single Chang-Refsdal microlens (e.g. Yonehara l2()()lll 
or a straight-line caustic (e.g. Shalvaoin 120011 Shalyapin 
et al. Eona. However, microlensing amplification maps 
(Wambsganss, Paczynski & Schneider llhtfOII even for com¬ 
ponent A, which is expected to have the smallest optical 
depth among the 4-images, suggest the possibility of a sub¬ 
stantial amount of caustic clustering effects. Therefore, a 
single microlens or almost straight-line fold caustic model 
may not be adequate to interpret the complex HAEs 
observed for the Einstein Gross. We have developed a 
consistent Local HAE Caustic Model (LOHCAM; Lee 
Eonai based upon N-body microlenses using a ray shoot¬ 
ing method (Kayser, Refsdal & Stabell I1986|l which si¬ 
multaneously takes into account the extended nature of 
the background source and the shear effects due to the 
macro deflector (e.g. lensing galaxy). LOHCAM implic¬ 
itly assumes that the HAEs are produced by those caus¬ 
tics located near the source and that the additional back¬ 
ground amplification due to the non-local microlenses may 
be considered constant and negligible during a local HAE; 
it makes use of the fact that the microlensing amplifica¬ 
tion in the rayshooting method depends only on the size 
of the source and on the strength of caustics, regardless of 
the unlensed flux of the source. 
































D.W.Lee et al: Possible evidence of an accretion disk with a central hole in QSO 2237+0305 


3 



Fig. 2. OGLE+GLITP data and model fitting results. Combined V-band light curves for image C and the possible 
caustic model, (a) The OGLE (231 points) & GLITP-ISIS (53 points; around day 1500) data are shown altogether 
with the simulated light curve (solid line) using a disk with a central hole source model, (b) The corresponding HAE 
caustic network (solid line) is shown with the straight track (dotted lines) for the moving source (Y=0). The concentric 
circles represent the sizes of the best-fitted inner and outer radii of the face on accretion disk. The sizes of the caustic 
network and of the source are represented in time units (days) along the X and Y axes. 


For source models, we restrict ourselves to consider 
only two simple models: uniformly bright circular disk 
with a hole and a uniformly bright disk models assum¬ 
ing face on status. According to the traditional unified 
theory of AGN, we can naturally imagine the presence of 
a central hole in an accretion disk (i.e. an accretion disk 
with a central hole; see Fig. 2 Thorne & Pricewith 
a SMBH in the heart of a quasar. The free parameters 
in LOHCAM are the mass and coordinates of the local 
microlenses, the inner and outer radii of a circular source 
assumed to be a face-on accretion disk with a hole sur¬ 
rounding a central SMBH and the background fiux level; 
i.e. 3N + 3 parameters. Additional pre-fixed parameters 
are the external shear, the shear direction angle defined 
with respect to the relative direction of the source motion 
and the number N of local microlenses. 

The unknown minimum background amplification is 
assumed to correspond to the minimum flux levels present 
in each of the A (see Fig. la around day 800) and C (see 
Fig. 2a around day 700) light curves. These levels are also 
found to be representative of minima recorded in the com¬ 
pilation of past data (Wyithe, Webster & Turner EOOOII . 
We adopt representative shear values of 0.42 and 0.62 at 
A & G from a general macro lens model for the case of 
the potential exponent (/3 = 1) (Witt, Mao & Schechter 
though the macrolensing parameters (e.g. ct, 7 ) 
are somewhat dependent on the adopted models (see also, 
Witt & Mao 1100^ Wambsganss & Paczvnski 1 199411 . We 


also note that although Schmidt, Webster & Lewis 
have reconstructed the macrolens model for the lensing 
galaxy of QS02237 by using the observed light distribu¬ 
tion (e.g. a bulge and bars), their model could not take 
properly into account the dark halo effect. Then, in order 
to take into account the dark halo of the lensing galaxy, 
we have chosen the macrolens model parameters from the 
most general model, the isothermal lens model (/3 = 1). 

Based on the normalized lens equation (Kayser, 
Refsdal & Stabell lTfiSfil Paczvnski ll986|l . 


V = 


1 + 7 0 

0 1-7^ ^ 


( 1 ) 


where 7 is the external shear, rrii represents the nor¬ 
malized mass of star i, 77 , ^ and the normalized co¬ 
ordinates in the source and the lens plane, respectively; 
the resulting microlensed lightcurves are then calculated. 
For simplicity, we ignore the totally unknown smooth 
dark matter term (ctc); this could slightly affect some 
of the scaling relations (Kayser, Refsdal & Stabell II 9861 
Paczvnski If 986(1 . Furthermore, due to the weak flux vari¬ 
ation seen in the D light curve (Wozniak et al. I2flflfl|l of 
QSO 2237+0305, we assume that the intrinsic flux varia¬ 
tions of the quasar are negligible and that all flux varia¬ 
tions are due to microlensing effects during the monitoring 
campaigns. 

When fitting the light curves, we tried to find the min¬ 
imum number N of local HAE microlenses in order to re¬ 
produce the whole light curves and to interpret the light 
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curves as smooth as possible where data are not avail¬ 
able. To convert the observed timescale into dimension¬ 
less normalized units for simulations, arbitrary time nor¬ 
malization is applied to the light curves, (e.g. 100 days). 
From the observed data, especially those of A, we notice 
that the average time scale between 2-consecutive turn¬ 
ing points in the light curve is about 100 days. This fact 
may imply an important clue for the sizes of the source 
and caustics. However, the exact value for the time nor¬ 
malization is not important, because it produces the same 
results only if the initial conditions are same in absolute 
scale (see Eq. 2). By introducing a time normalization, 
the observed light curves transform into normalized (di¬ 
mensionless) scale light curves. Then we can make use 
of the dimensionless normalized microlensing simulations 
compared with the normalized light curves without con¬ 
sidering any cosmological effects. 

Of course, when we convert a normalized unit into a 
real physical unit, then we should take into account cosmo¬ 
logical effects. A large number of preliminary simulations 
were then performed to identify plausible ranges for the 
initial parameters. Starting from the binary lens plus ex¬ 
ternal shear cases ILee lHfflTl Lee, Chang & Kim lHhk^ . we 
have searched various caustic cuts to match the observed 
light curves. However, we could not find any satisfactory 
results from the binary lens model, so we increased the 
number of microlenses to find the possible minimum num¬ 
ber caustic solutions compatible with the observational 
data. Finally, we found that the minimum number of mi¬ 
crolenses needed to account for the observed HAEs is 5 
for both images A & C (Figs, lb & 2b). 

Next, applying a simultaneous minimization method, 
the downhill simplex algorithm (Press et al. II !TO2II . the 
best chi-squared set of parameters characterizing the ob¬ 
servational data was finally derived. It is important to 
note that knowledge of the apparent source velocity with 
respect to the caustics directly determines the absolute 
sizes of the source and caustics and, consequently, the 
mass of each microlens through the time normalization 
of LOHCAM. Through the LOHCAM analysis, we can 
directly derive the physical mass (M) of a microlens based 
on the normalized mass as follows; 

Dt 

M = — — X (nm X lOOdays x Vt)"^ (2) 

4G DsDls 

where c is the speed of light, G is the gravitational con¬ 
stant, Dl, Dls^Ds represent angular diameter distances, 
nm represents the normalized mass, lOOdays represents 
the basic normalization unit in our simulations, Vr is the 
projected transverse velocity of the source with respect to 
the caustics in the source plane measured by an observer. 

4. Results 

Although it is extremely hard to prove the uniqueness 
of our proposed caustic models as we never observe the 



Fig. 3. The difference between two source models for the 
A light curve. The upper graph shows the case optimized 
for a disk with a hole source model (y^ = 1.72). The lower 
graph shows the case optimized for a disk source model 
(slightly offset for clarity; = 2.38). Note the much bet¬ 
ter matches for the case of the disk with a hole model for 
the double peak, turning points and a small dip. These 
two curves are the best-optimized cases by using a least 
number of microlenses (i.e. 5-microlenses), respectively. 


caustics, we could derive some interesting and consistent 
results for the microlensing effects in the Einstein Cross. 
From the numerical simulations, we find that only a source 
model with a hole inside may easily account for the dou¬ 
ble peak observed in the light curve of component A (Figs, 
la & b). Without such a hole, the simulated light curves 
fail to reproduce this double peak (Fig. 3). Even if we in¬ 
clude additional microlenses (6 or 7 in total), we failed to 
find any good fitting results for the double peak with a 
disk source without a hole. We note that this double peak 
modeling for A is very similar to Grieger et al. dUHi ’s the¬ 
oretical prediction for the case of an accretion disk with a 
central hole source. This result thus provides the first pos¬ 
sible evidence of an accretion disk with a central hole for 
the far-UV continuum source of a distant quasar, namely 
QSO 2237+0305, observed at a redshift z=1.695 in the 
V band. In addition, our modeling produces a remark¬ 
ably good fitting over the entire data points (Figs, la & 
2a), not just over selected parts of HAEs ('Youeha,ra l2001l 
Shalyapin 120011 Shalyapin et al. 12003). Our results are 
directly understood from the constant motion of the face 
on accretion disk source along a straight track across the 
caustic networks (Figs. 1 & 2). The final fitting val¬ 
ues are small, but not sufficiently small. This may be in¬ 
duced by our simple uniform brightness modeling for the 
surface brightness distribution of the source. As we can 
see from the fitting results, most of parts of the A data 
are well reproduced by the LOHCAM. Though the double 
peak modeling is evident, at the double peak region there 
is a slight misalignment between our modeling and the 
data; we need to make more sharp peaks. If some one can 
take into account more realistic surface brightness distri- 
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Table 1. The best optimized model fitting results through LOHCAM. 



A 

C 

Inner source radius ((Vt/ 25, 000fcms~^)Af/) 
Outer source radius ((Vt/ 25, 000fcms“^)A?7) 
Masses (10“^ x (Vt/ 25, OOOfcms“^)^M0) 
Reduced x^ 

254 (+150/-119) 

1318 (+23/-73) 

1.07, 1.91, 1.19, 1.23, 1.29 
1.72 

< 349^ 

1297 (+144/-77) 

102.02, 2.04, 1.20, 0.75, 0.34 
1.86 


The inner source radius for the case of C is not well constrained. This may be due to the lack of data around the peak 
and the fact that the peak of HAE of C image is already affected by another large caustic (see Fig. 2). Quoted errors for 
the source sizes are la. The errors were calculated from the 2-D (i.e. inner source radius versus outer source radius) Ax^ 
contour map corresponding to Icr level, while fixing the other parameters at their best-fitted values. Please note that these 
values are functions of the transverse velocity and these values are derived assuming Vt = 25, OOOfcms”^ in the source plane. 


bution of the source with our caustic and source model, 
the corresponding results should be near to unity be¬ 
cause most of loss occurs at the peak region in our 
modeling. Shalyapin et al. II2()()2II made a HAE fitting for 
the light curve of A only based on GLITP-PSF data and 
derived a small x^- However, we note that GLITP-PSF 
data has larger error bars compared to GLITP-ISIS data; 
so it is easy to have a smaller x^ value. Moreover, unlike 
any other HAE studies made so far, our model fitting pro¬ 
duces independent and similar solutions for the source size 
from both the A & G light curves (Table 1). In our opinion, 
this testifies the coherence of the LOHGAM results, also 
suggesting that the contribution of the unknown smooth 
dark matter is not sufficient to produce over-focused cases 
in this lens system (Kayser, Refsdal & Stabell ITUHtill . 

From the best fitting results, we may now infer the 
physical sizes for the accretion disk and the central hole of 
the quasar and the masses of the microlenses. Please notice 
that all the results of LOHGAM only depend on the trans¬ 
verse velocity of the source with respect of caustics in the 
source plane measured by an observer; the size of source 
is linearly proportional to the transverse velocity and the 
masses of microlenses are proportional to the square of the 
transverse velocity. Although Kayser & Refsdal lllbbbll es¬ 
timated 6, 000 A:tos“^ for the projected transverse velocity 
of the source with respect to the caustics in the source 
plane, recently Kochanek li™il calculated a more real¬ 
istic transverse velocity of QS02237 in the source plane 
based on statistical analysis of the OGLE light curves. His 
estimation is roughly between 10,000 ~ 40, OOO/cms”^. 
Adopting a mean value (Vr = 25, 000kms~^) of his esti¬ 
mation, we estimate the size of source and the masses of 
microlenses (Table 1). 

For the cosmological model, we assume No = 
65kms~^Mpc~^ for the Hubble parameter and a flat 
Universe with a matter density flm = 0.3. The best- 
fitted outer and inner radii of the source for A are then 
found to be 1318(+23/ - 73) x (Vt/25, 000A:tos- 1)A{7 
and 254(+150/ - 119) x (Ut/25, 000A:tos- 1)AU, respec¬ 
tively. Similar values are obtained from the independent 
analysis of C (Table 1). From the derived size of the in¬ 
ner radius (Rm) of the accretion disk and using a rela¬ 
tion for the last stable orbit {Rin ~ ^Rsw) of a black 


hole with a Schwarzschild radius (Rsw), the mass of the 
SMBH in QSO 2237+0305 is estimated to be in the 
Icr range 4.3(+2.5/ - 2.3) x (Ut/25, OOOfcms-QlO^M©. 
This mass estimate is the largest one compared to known 
SMBHs (Gebhardt et al. 1200011 . Our mass estimation 
for the SMBH of QS02237 is almost compatible with 
Kochanek’s recent result and 10 or 100 times 

higher than Goicoechea’s result II2003II . if we adopt Vr = 
25,000fcms“^. 

When adopting Vt = 25, OOO/cms”^, the estimated 
mass range for the microlenses is of the order of lO^^M©, 
although one star with a mass IM© is required to fit the 
HAE of G (see Table 1). Then, our mass estimation of mi¬ 
crolenses in the Einstein Cross is almost consistent with 
the results by Kochanek But, if we adopt a smaller 

transverse velocity, Vr = 6 , 000A:tos“^ (Kayser, Refsdal & 
Stabell lTOSbl Kayser & Refsdal [198911 . the estimated mass 
range is between 10“"'^M© ~ 10“^M©. This is below the 
hydrogen-burning limit; it can thus represent a range of 
masses compatible with extra-galactic brown dwarves or 
free floating Jupiter-like planets. Refsdal & Kayser II1993II 
have previously anticipated that microlensing by small 
masses with a large source size {Re < Rs) could explain 
the observed variability in the Einstein Cross. Our results 
may support their claim on the point that the source size 
is larger than the Einstein length of a mean mass of mi- 
crolenses, since our source size estimation is slightly larger 
than the Einstein length of mean microlenses (Fig.l & Fig. 
2). It also suggests the possibility that there can be numer¬ 
ous planetary masses in the distant lensing galaxy, if we 
adopt Vt = 6, 000kms~^. The exact masses of microlenses 
in the Einstein Cross still remain uncertain since there are 
very different estimations for the transverse velocity be¬ 
tween 6, OOOfcms”^ ^ 40,000fcms“^. Accordingly, for ex¬ 
ample, the mass range corresponding to 100 days in the 
source plane is estimated between 0.0007A/© ^ 0.033M© 
depending on the adopted transverse velocity. Though the 
exact masses of microlenses and the size of source are func¬ 
tions of the transverse velocity, the derived caustic size 
and the source size are invariant in time dimension or in 
normalized unit in LOHCAM (Fig. 1 & Fig. 2). 

Although we have firstly estimated the masses of mi¬ 
crolenses in this galaxy under the assumption of the trans- 
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verse velocity, it is not possible to derive or compare the 
expected surface mass density from our model. Because 
we must know both the mass and the spatial density of 
microlenses in order to derive the surface mass density, 
we can not estimate the surface mass density as we can 
not know the mean spatial density of microlenses from our 
modeling; we only estimate a few masses of microlenses, 
not their mean spatial density. Furthermore, the surface 
mass density in such a locally small region (see Fig. 1 & 
Fig. 2) can be very fluctuated from the value expected 
in macro-lens models. However, it seems that the spatial 
density of microlenses in this system is not quite dense, 
otherwise our simple 5-microlens model could not repro¬ 
duce the observed HAE variations. 

Our current estimate for the UV-continuum source size 
of the quasar measured in V-band (e.g diameter ~ 3.94 x 
10 ^®(Vt/ 25, 000fcTOS“^)cTO) is slightly larger than the es¬ 
timated values by Shalyapin et al. COP and Yonehara 
lj2()()l|l . if we adopt Vt = 25, 000A:tos“^. Although Wyithe 
et al. Il2()()()ll argued that the continuum source size of 
QSO 2237 would be much smaller than the Einstein length 
for the mean average mass of microlenses and neglected 
the low mass hypothesis by Refsdal & Stabell , our 

results clearly show that it is possible to generate the 
strongest HAEs observed by OGLE with the combina¬ 
tion of small size caustics compared to the source size. 
We also note that Wyithe et al.’s II2()()()II prediction for 
the HAE in the Einstein Cross is completely failed; they 
predicted too much strong HAEs and it seems that they 
have based on too small source or too large caustics. If 
we assume a small transverse velocity, Vt = 6 , 000A:tos“^, 
they should be planetary masses. Then, our results are 
supporting Refsdal & Kaiser ’s idea that there can 

be a large amount of planetary masses, though the exact 
masses are still pending on the true transverse velocity in 
the source plane. Therefore, studying the transverse ve¬ 
locity in the source plane with a statistical analysis for 
the future OGLE light curves is very important so as to 
determine the exact size of the source and masses of mi¬ 
crolenses in this system. 

Considering two mutually perpendicular source track 
directions with respect to the galaxy center, we find that 
the whole HAEs of A are better modeled with the direc¬ 
tion parallel to the shear direction (i.e. defined toward 
the galaxy center), whereas the whole HAEs of C are 
smoothly fitted with a track direction perpendicular to 
the local shear (Figs. 1 & 2). Given the orthogonal posi¬ 
tions of A & C with respect to the galaxy center, these 
results may provide a rough, but consistent clue for the 
direction of the relative transverse motion of QSO 2237. 
Assuming that the random motions of the individual mi¬ 
crolenses are negligible, we cautiously conclude that the 
direction of transverse motion of the source is more likely 
roughly parallel to the line connecting A & B, instead of 
C & D. This finding may imply that the true transverse 
velocity of the quasar can be larger than expected since 
the predicted observer’s motion is parallel to the C-D di¬ 
rection (Witt & Mao ll99lll . 


Assuming that our interpretation of the source motion 
in the sky is correct, we can predict that such abrupt small 
fluctuations (cf. around day 1600 in Fig. 2) will be more 
frequent in image G, whereas image A will show smoother 
variations when microlensing HAE occurs. Also, the peak 
of HAE will be more sharp in C light curves as the pre¬ 
dicted source’s moving direction is perpendicular to the 
shear direction. Therefore, our current interpretation of 
the source motion should be further investigated on the 
basis of new data for the Einstein Gross. Especially, mon¬ 
itoring observation of the centroid shift in QS02237 by 
the Space Interferometry Mission (SIM) is of great inter¬ 
est for us. If small mass lenses are dominant in this system, 
we could not expect to observe a significant centroid shift 
even during evident HAEs down to micro arcsecond scales. 
So the future SIM project may provide valuable informa¬ 
tion about the true transverse velocity and the masses of 
microlenses in this system. 

Our new technique, the LOHCAM, will play a key role 
for the future HAE analysis for extra-galactic microlensing 
events. 
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